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The present knowledge of Lamb shift, fine- and hyperfine structure of the 2S and 2P states in 
muonic deuterium is reviewed in anticipation of the results of a first measurement of several 2S — 2P 
transition frequencies in muonic deuterium (/id). A term-by-term comparison of all available sources 
reveals reliable values and uncertainties of the QED and nuclear structure-dependent contributions 
to the Lamb shift, which are essential for a determination of the deuteron rms charge radius from 
/rd. Apparent discrepancies between different sources are resolved, in particular for the difficult 
two-photon exchange contributions. Problematic single-sourced terms are identified which require 
independent recalculation. 


I. INTRODUCTION 

Laser spectroscopy of 2S —i 2P Lamb shift transitions 
in muonic atoms and ions promises a tenfold improve¬ 
ment in our knowledge of charge and magnetic radii of 
the lightest nuclei (Z = 1,2 and higher). Our recent mea¬ 
surement mm of the 2S Lamb shift and the 2S hyperfine 
splitting (HFS) in muonic hydrogen, /rp, in combination 
with accurate theoretical calculations by many authors, 
summarized in Ref. [3], has revealed a proton root-mean- 
square (rms) charge radius of 

r p = 0.84087 (26) exp (29) theo fm = 0.84087 (39) fm. 

(1) 

This is an order of magnitude more accurate than 
the value of r p = 0.8775(51) fm evaluated in the CO¬ 
DATA least-squares adjustment [4] of elastic electron- 
proton scattering UM and many precision measurements 
in electronic hydrogen [7j. 

Most strikingly, however, the two values differ by 7 
combined standard deviations (7cr). Despite numerous 
attempts in recent years to explain this “proton radius 
puzzle”, it remains a mystery mm- Taken at face value, 
this 7 a discrepancy constitutes one of the biggest discrep¬ 
ancies in the Standard Model. Further data are clearly 
required to shed light on this puzzle. 

Muonic deuterium, fid, has been measured in the same 
beam time as /ip mu, and the data are now nearing pub¬ 
lication [1C)|. We anticipate here that the experimental 
accuracy of the various 2S — 2P transitions is of the or¬ 


der of 1GHz, or, equivalently, ~ 0.004meV[(] Ideally, 
theory should be accurate on the level of 0.001 meV to 
exploit the experimental precision, and to determine the 
deuteron charge radius, r dl with tenfold better accuracy, 
compared to the CODATA value [3] 

r d (CODATA) = 2.1424(21) fm. (2) 

The CODATA value originates from a least-squares ad¬ 
justment of a huge amount of input values, such as 
the deuteron charge radius from elastic electron scatter¬ 
ing mm 

r d (e — d scatt.) = 2.130(10) fm, (3) 

but also the proton radius from electron scattering mm. 
These radii are connected because the CODATA adjust¬ 
ment includes many transition frequencies in hydrogen 
(H) and deuterium (D) [31 tZ]- In particular, the squared 
deuteron-proton charge radius difference, 

r 2 - r 2 = 3.82007 (65) fm 2 (4) 

is known with high precision from laser spectroscopy of 
the isotope shift of the IS — 2S transition in electronic 
hydrogen and deuterium }l3j , and state-of-the-art the¬ 
ory [2]. Using Eq. © and the muonic hydrogen proton 
radius given in Eq. (|T|) we determined a value of j2l 

r d (muonic r p ) = 2.12771 (22) fm. (5) 

Note that the discrepancy of the deuteron charge radii 
given in Eq. ([2| and Eq. ([5]) is not a new discrepancy, 
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but rather a result of the proton radius discrepancy: Both 
values of the deuteron radius depend on the isotope shift 
in Eq. Q. Hence, discrepant values of the proton radius 
will result in discrepant values of the deuteron radius. 

The upcoming /id data m, on the other hand, will 
provide a “muonic” value of the deuteron radius that is 
independent of the proton charge radius. As such, it will 
shed new light on the proton radius puzzle. 


We anticipate here that the theory of the Lamb shift 
in muonic deuterium is limiting the accuracy of the 
deuteron charge radius from /id, mainly due to the un¬ 
certainty of the deuteron polarizability contribution of 
0.020 meV which corresponds to a relative uncertainty 
of 1%. Nevertheless, the deuteron charge radius from 
/rd m will have a nearly three times smaller uncertainty 
than the current CODATA value (Eq. <§)■ 

To put this uncertainty of 0.020 meV into another per¬ 
spective: The “proton radius puzzle” in muonic hydro¬ 
gen, when expressed as a “missing part” in the theory of 
muonic hydrogen, amounts to 0.329 meV. 


This article is organized as follows: We first sum¬ 
marize the current knowledge of the muonic deuterium 
Lamb shift theory (Sec. Ill) which is required to deter¬ 
mine the deuteron rms charge radius r a from the pd 
measurement m • We separate the Lamb shift theory 
into “radius-independent” terms that do not depend on 
the nuclear structure (Sec. Ill A), terms that depend ex¬ 
plicitly on the rms charge radius (Sec. Ill B), and the 
deuteron polarizability contribution that constitutes the 
main theoretical limitation (Sec. Ill C). In Sec. |IV| we 
list all contributions to the 2S hyperfine splitting (HFS) 
in pd. The 2S HFS depends on the magnetic properties 
of the deuteron through the Zemach radius. Other nu¬ 
clear structure contributions matter, too, so we separate 
again terms: Sec. |IVB] lists the terms that do not depend 
strongly on the deuteron structure, Sec. |IV C'| is devoted 
to the Zemach correction, Sec. IV D is concerned with the 
deuteron polarizability contribution which has recently 
been calculated for the first time m This term con¬ 
stitutes the main uncertainty for the 2S HFS prediction. 
Additional contributions to the 2S HFS are mentioned in 
Sec. |IV E] The 2P fine structure is summarized in Sec. [7] 
and the 2P fine- and hyperfine level structure, including 
level mixing, is given in Sec. |VI| 


The sign convention in this article is such that the fi¬ 
nal, measured energy difference (Lamb shift, 2S-HFS, fine 
structure) is always positive. For the fine and hyperfine 
splittings this convention is the natural choice adopted 


by all other authors, too. For the Lamb shift, however, 
some authors calculate 2S level shifts and their published 
values have the opposite sign. This is because the 2S 
level is lower (more bound) than the 2P level (due to the 
dominant vacuum-polarization term item ffl in Tab. [jj, 
see Fig. [I] and positive level shifts decrease the measured 
2P-2S energy difference. The numbers we quote are all 
matched to our sign convention. 

Item numbers ff in the first column of Tab. |T] and 
Tab. |v| follow the nomenclature in Ref. [3|. In the ta¬ 
bles, we usually identify the “source” of all values en¬ 
tering “our choice” by the first name of the (group of) 
authors given in adjacent columns (e.g. “B” for Borie). 
We denote as average “avg.” in the tables the center of 
the band covered by all values under consideration, 
with an uncertainty of half the spread, i.e. 

avg. = ^ [MAX(uj) +MIN(uj)] 

± g [ MAX (^) ^ MIN( Wi )] 

Throughout the paper, Z denotes the nuclear charge 
with Z = 1 for the deuteron, a is the fine structure con¬ 
stant, m r is the reduced mass of the muon-deuteron sys¬ 
tem. “VP” is short for “vacuum polarization”, “SE” is 
“self-energy”, “RC” is “recoil correction”. “Perturbation 
theory” is abbreviated as “PT”, and SOPT and TOPT 
denote 2 nd and 3 rd order PT, respectively. 


II. OVERVIEW 


The n = 2 levels in muonic deuterium are sketched in 
Fig. [3 The Lamb shift, i.e. the splitting between the 2S 
and the 2P!/ 2 state, is sensitive to the rms charge ra¬ 
dius of the deuteron, as detailed in Sec. |III| In contrast, 
2S hyperfine splitting (HFS) is caused by the magnetic 
interaction between the muon spin and the magnetic mo¬ 
ment of the deuteron. The finite deuteron size results in 
a finite magnetization distribution inside the deuteron, 
and makes the 2S HFS depend on the so-called Zemach 


radius of the deuteron, as explained in Sec. IV 


The first calculation of the Lamb shift in muonic deu¬ 
terium was published by Carboni m in 1973. More 
elaborate calculations of QED effects in muonic atoms 
were introduced with the seminal paper by Borie and 
Rinker flTj in 1982. 

Later, Pachucki [TS] and Borie [TS] presented very de¬ 
tailed calculations of many terms for muonic hydrogen. 
Borie then extended her pp calculations [TS] to the case 
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FIG. 1. The 2S and 2P energy levels in muonic deuterium. The inset on the right displays the shifts A 1/2 and A 3/2 of the 


2P (F = 1/2) and 2P(.F = 3/2) levels due to the mixing of the F = 1/2 and F = 3/2 states, respectively, as described in Sec. VI 
The figure is not to scale. 


of muonic deuterium |20j . After our measurements in 
muonic hydrogen m and deuterium, Borie revisited the 
theory of the n = 2 energy levels in light muonic atoms 
(/rp, /id, /i 3 He, and /r 4 He) in Ref. [ 21 ] • The published pa¬ 
per [21, (available online 6 Dec. 2011) has subsequently 
been superseded by the arXiv version of the paper [22] . 
At the time of this writing, Borie’s paper on the arXiv 
has reached version 7 (dated 21 Aug. 2014). This is the 
first source of knowledge on /id summarized in here. 

The second source is the group around Faustov, Kru¬ 
tov, Martynenko etal ., termed “Martynenko” in here for 
simplicity. They have published an impressive set of pa¬ 
pers on theory of energy levels in light muonic atoms. 
At the time of this writing, the 2011 paper [231 was the 
most recent one on the Lamb shift in /id, and we based 
our summary on this paper. Later, Ref. [24] from the 
Martynenko group appeared, with only minor differences 
in the results compared to Ref. [23] . For simplicity, we 
still base our compilation of Lamb shift contributions on 
the earlier, more detailed, paper [23]. In particular, equa¬ 
tion numbers and table entries refer to Ref. [23], unless 
otherwise noted. For the 2S HFS, Ref. m is the main 
source of numbers from the Martynenko group. 

After the advent of the proton radius puzzle from 
muonic hydrogen, many groups have revisited and im¬ 


proved the theory on muonic hydrogen in an (unsuc¬ 
cessful) attempt to identify wrong or missing theory 
terms large enough to solve the puzzle (see our com¬ 
pilation [3 for a detailed overview). Thankfully, two 
groups have (re-)calculated many terms not only for the 
case of muonic hydrogen, but also for muonic deuterium 
(and /UHe and /x 4 He): Jentschura, and Karshenboim’s 
group with Ivanov, Korzinin, and Shelyuto, have pub¬ 
lished many papers on muonic deuterium which are in¬ 
cluded in the present compilation. 

III. LAMB SHIFT IN MUONIC DEUTERIUM 
A. QED contributions 

One-loop electron vacuum polarization (eVP) (Fig.j2|, 
the so-called Ueliling term [25], accounts for 99.5 % of the 
nuclear-structure-independent part of the Lamb shift in 
/rd. It is therefore mandatory to double-check this term 
as thoroughly as possible. 

Borie has argued H3 mm that the Uehling term 
should ideally not be treated perturbatively. Instead, the 
Dirac equation should be solved numerically after adding 
the Uehling potential to the electrostatic Coulomb po¬ 
tential. For light muonic atoms such as muonic deu- 
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FIG. 2. Item #1, the leading order 1-loop electron vacuum 
polarization (eVP), also called Uehling term. 

terium, however, both approaches should give accurate 
results |22| . This has been demonstrated for muonic 
hydrogen , where the nonperturbative result of Indeli¬ 
cate |26j is in excellent agreement with the perturbative 
results of Pachucki ubi and Borie KHi, see Ref. [3]. 

For muonic deuterium , only perturbative calculations 
exist, albeit with two slightly different approaches: 

Martynenko etal. [23], Jentschura E3 ng and 
Karshenboim etal. :29 ;i [30] calculate the leading order 
eVP nonrelativistically (item #1 in Tab. [I]), and apply 
a relativistic correction (item #2). The most important 
item #1 is in excellent agreement for all three authors, as 
well as with the value of 227.635 meV obtained by Car- 
boni m in 1973. For item #2 see below. 

Borie, in contrast, uses relativistic Dirac wave func¬ 
tions to evaluate the relativistic Uehling term (item #3). 
The relativistic recoil correction to eVP of order a(Za) 4 
(item #19) has to be added, to be able to compare all 
four results. 

It is very reassuring that these results are in excellent 
agreement, with one exception: Item #2, rel. corr. 
(Breit-Pauli), from Martynenko [35] . 0.0177meV, dif¬ 
fers from the value of 0.02178 meV, calculated by the 
other three groups, who agree: Borie [22] Tab. 1, 
Jentschura [3Tj Tab. I and [28] Eq. (17), and Karshen¬ 
boim [30] Tab. IV. This item #2 should be the 
sum of Martynenko’s rows 7 and 10 ( relativistic and 

(a) (b) (c) 





FIG. 3. Item #4, the two-loop eVP (Kallen-Sabry) contribu¬ 
tion. 


FIG. 4. Item #5, the one-loop eVP in 2-Coulomb lines. 

VP corrections of order a(Za) 2 in first and second 
order PT). For muonic helium-3 and -4 ions [ 32 ] 
the sum agrees exactly with the numbers given by 
Jentschura in Ref. m Eq. (17) and by Karshenboim in 
Ref. [ 30 ] Tab. IV. Martynenko confirmed that their value 
for our item #2 for muonic deuterium (their rows 7 and 
10 in Tab. I of Ref. [ 23 ]) contains an error [ 331 . 

The average of items #l+#2 or #3+#19 is thus calcu¬ 
lated from the other three sources only, with excellent 
agreement: 

AE (one — loop eVP with rel. corr.) 

= 227.65658 ± 0.00020 meV (7) 

Our item #4, the two-loop electron-VP correction , usu¬ 
ally called “Kallen-Sabry” contribution [33], displayed in 
Fig- i is the second largest “purely QED” contribution 
to the Lamb shift. Borie [22] and Martynenko [23] give 
values which are in very good agreement. 

Our item #5, the one-loop eVP insertion in 2 Coulomb 
lines shown in Fig. [4] has been calculated with very 
good agreement by Borie [ 22 ], Martynenko [ 23 ], and 
Jentschura |28J. 

Karshenboim et al. :29] give the sum of items #4 and 
#5. It is in agreement with the sums from Borie and 
Martynenko. We use Karshenboim’s value because it is 
given with more significant digits. 

(a) (b) (c) 




FIG. 5. The three contributions to Light-by-light scatter¬ 
ing: (a) Wichmann-Kroll or “1:3” term, item #9, (b) Vir¬ 
tual Delbriick or “2:2” term, item #10, and (c) inverted 
Wichmann-Kroll or “3:1” term, item #9ah 

























5 


Light-by-light (LbL) scattering (see Fig. [5]) contains 3 
terms, Wichmann-Kroll, or “1:3” LbL (item #9 in [3]), 
Virtual Delbriick, or “2:2” LbL (item #10 in [3]), and the 
“inverted Wichmann-Kroll” or “3:1” LbL (called “new” 
item in 0 ). For definiteness, we label this term “item 
#9a” from now on. Considerable cancellations occur in 
the sum of all three terms which has been evaluated in 
Ref. [35]. Both Borie and Martynenko calculate only #9, 
and adopt the full result from Karshenboim [35]. We use 
Karshenboim’s result [35] . 

Item #11 muon self-energy (SE) correction to elec¬ 
tron vacuum polarization (eVP) a 2 3 (Za) 4 is displayed 
in Fig. § Jentschura [28] (Eq. (29b), Fig. 2) and 
Karshenboim [37] (Tab.VIII (a), Fig. 6a) agree in the re¬ 
sult for the complete calculation, -0.00306 meV. Marty¬ 
nenko [23] Eq. (80) calculates only the contribution from 
Fig. [6](a) , following Pachucki’s Eq. (39) in Ref. [18] • Also 
Borie [22] calculates part of this term in her Appendix C. 

Higher order corrections to the muon self-energy and 
vacuum polarization are denoted items #12, #13, #21, 
#30* and #31* in Tab. [TJ For muonic hydrogen t 3j we 
used Borie’s value of item #21, noting that this includes 
item #12. Afterwards, Karshenboim etal. m have re¬ 
calculated many of these small terms. We construct the 
corresponding sum from each source, which we average. 

Item #12 is shown in Fig. [7] It is Martynenko’s item 
29. This contribution has been confirmed by Karshen¬ 
boim [37] Tab. VIII (d). As mentioned in Ref. 0], 
item #12 is included in Borie’s value for item #21. 

Item #13, mixed muon-electron VP is depicted in 
Fig. i Borie and Martynenko calculated only the 
contribution from Fig. |8](a) , see Fig. 3 in Ref. [38]. This 
is Martynenko’s item 3, “VP and MVP contribution in 
one-photon interaction”. Karshenboim gives the sum of 
both diagrams in Fig. [8] in Ref. [37] Tab. VIII (d). 


(b) 


(c) 






FIG. 6. Item #11, muon self-energy corrections to the elec¬ 
tron vacuum polarization a 2 (Za) 4 . This figure is Fig. 2 from 
Jentschura [36]. It corresponds to Fig. 6(a) from Karshen¬ 
boim /57. 


Item #30* (#31*) is somewhat similar to item #12 
(#13), with the electron (muon) loop replaced by a 
hadronic VP loop, see Fig. Q. It has only been cal¬ 
culated by Karshenboim [37], Tab. VIII item (e) ((c)). 

Item #32, the muon VP in SE correction shown in 
in Fig. 0 is not included as a separate item in our 
Tab. [T] It should already be automatically included in 
any QED value which has been rescaled from the QED 
of electronic deuterium by a simple mass replacement 
m e —1 [35] • The size of this item #32 can be esti¬ 

mated from the relationship found by Borie (HI], that 
the ratio of hadronic to muonic VP is 0.66. With 
Karshenboim’s value of item #30 * m one would get 
A£’(#32) = —0.000024/0.66 meV = -0.000036 meV. 

Item #21, higher order correction to pSE and pVP, is 
Borie’s muon Lamb shift, higher orders, calculated in her 
Appendix C of Ref. [22]. This item includes item #12, 
eVP loop in self-energy a 2 (Za) 4 , as explained on p. 131 
of Ref. 0 0 

The sum of items #12, #13, #21, #30*[^J and 
#31* agree well enough to justify taking the average, 
—0.00178 ± 0.00014 meV as our choice. 

Item #14, hadronic VP, is evaluated by Borie [22] 
(p. 5) as 0.013 meV, who assigns a 5% uncertainty to 
this estimate which is based on Refs. [mgDHig. 
Martynenko’s value of 0.0129 meV [23], Tab. 1, row 31, 
agrees very nicely. They quote Ref. [44] 03], and estimate 
the uncertainty to 5% as well. 

The previous items #15 and #16 in Ref. [3] are higher 
order corrections to the hadronic VP and have only 
been calculated for muonic hydrogen by Martynenko’s 



FIG. 7. Item #12, eVP loop in SE are radiative corrections 
with VP effects. This is Fig. 11 from Martynenko [23] which 
is the same as Fig. 4 in Pachucki m- It is Karshenboim’s 
Fig. 6(d) in Ref. [37]. 


2 There is a typo in footnote f of Tab. 1 in Ref. [3], where we 
wrote, item #12 “is part of #22” (instead of #21). 

3 The asterisk * indicates that this item had not been considered 

for muonic hydrogen in Ref. [3]. 
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(a) (b) 




FIG. 8. Item #13, the mixed eVP-/zVP contribution. 

group [44l [45] , where they are small (0.000047 meV and 
-0.000015 meV, respectively). It is expected that their 
magnitude will be similar in muonic deuterium. Hence 
these terms are included in the 5% uncertainty assigned 
to item #14. 

Item #18 in Ref. [3] is the recoil finite size contribution 
in Borie EH 121]. According to Pachucki, this item #18, 
which was first calculated by Friar [46], should be dis¬ 
carded 07, as we did for muonic hydrogen in Ref. (3]. 

Item #17 is the Barker-Glover correction [3S], termed 
additional recoilby Borie. This item includes the Darwin- 
Foldy (DF) term that arises from the Zitterbewegung of 
the nucleus. For a spin-1 nucleus such as the deuteron 
(as well as for the spin-0 4 He nucleus) this DF term is 
absent [49]. Different conventions are used in the liter¬ 
ature [5f)U52| which has caused some confusion (see Ap¬ 
pendix |A|. As in the case of muonic hydrogen [3] (where 
the DF term is nonzero), we follow the “atomic physics” 
convention [52], also adopted by CODATA-2010 [4]. In 
this way, the charge radii from muonic hydrogen mm and 
deuterium m are directly comparable to the CODATA- 
2010 values (4], as well as the electronic H-D-isotope shift 
given in Eq. Q, all of which follow the same conven¬ 
tion miii]. 

Item #28^ is the rad. (only eVP) RC a(Za) 5 labeled 
“new” in Ref. 3]. For definiteness, we enumerate it as 
item #28^. It is the sum of three individual parts which 
sum up to 0.000093 meV in Ref. [21], Eq. (46b). 
Martynenko’s row 26 of Tab. I in Ref. [23], recoil corr. 



d 

FIG. 9. Item #30*, hadronic VP in SE contribution, corre¬ 
sponds to Fig. 6(e) in Karshenboim | I37J . 


(a) (b) 




FIG. 10. Item #31*, the mixed eVP- and hadronic VP con¬ 
tribution, comes from the Uehling correction to the hadronic 
VP correction. See Fig. 6(c) in Karshenboim 1371 . 

to VP of order a(Za) 5 (seagull term) is only the seagull 
term from the three terms evaluated by Jentschura, taken 
from Ref. [21], Eq. (29). We take the full result from 
Ref. [22], Eq. (46b). 

Item #24 are radiative recoil corrections of order 
a{Zaf and ( Z 2 a)(Za ) 4 , first introduced for /ip by 
Pachucki m, Eq. (51), based on Ref. [S3]. Borie writes 
(p. 9 of Ref. (21) that these terms correspond to Tab. 9 
and two additional terms from Tab. 8 of the review of 
Eides etal. [42j . The sum is -0.00302 meV for /zd. Mar¬ 
tynenko etal. (Ref. El, row 27 of Tab. 1, and Eq. (71)), 
on the other hand, evaluate only the terms from Tab. 9 
of Eides El, which gives -0.0026meV. We use Borie’s 
complete result. 

Item #29*, the a 2 (Za) 4 m contribution to the Lamb 
shift, is new and has not been considered for muonic hy¬ 
drogen in Ref. [1- We keep enumerating the items, mak¬ 
ing this item #29* [f] 

Martynenko gives this as the sum of rows 8 and 11 in 
Tab. 1 of [251 . relativistic and two-loop VP corrections of 
order a 2 (Za) 4 in 1st and 2nd order PT. They sum up to 
—0.0002 + 0.0004 = 0.0002 meV. Karshenboim calculated 
the complete correction of this order, with recoil correc¬ 
tions included and calls it eVP2 in Tab. VIII of Ref. [37] , 
Their value of 0.000203 meV replaces Martynenko’s rows 
8 and 11 [39]. 



FIG. 11. Item #32, muon VP in SE contribution, is automat¬ 
ically included in a rescaled electronic deuterium QED value 
of higher order SE contributions (see text). 
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Numerically this item #29* is of little practical impor¬ 
tance, because it is so tiny: 0.000173 meV in /up, and 
0.000203 meV in /ud. Of course, the calculation of this 
term was an important confirmation that previously un¬ 
calculated higher order terms are not responsible for the 
proton radius discrepancy. It is very reassuring that 
the two different approaches of Martynenko [23| and 
Karshenboim m give the same result. Interestingly, 
there is no such an agreement for the cases of muonic 
helium-3 and -4 ions. In view of our recent measure¬ 
ments in muonic helium-3 and -4 |54j , this disagreement 
may deserve further study, even though the size of the 
terms is small compared to the overall uncertainty. 

The sum of all contributions without explicit nuclear 
structure dependence summarized in Tab. |T| amounts to 

A£ r L a S d._indep. = 228.77356 ± 0.00075 meV. (8) 
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B. Radius-dependent contributions to the Lamb 
shift 


The radius-dependent contributions to the Lamb 
shift [22EH.3H are listed in Tab. Ill Generally, the finite 


size of the nucleus affects mainly the S states whose wave 
function is nonzero at the origin, where the nucleus re¬ 
sides. The main finite size contributions to the nS states 
have been given to order ( Za ) 6 by Friar |l(ij 


Ai£fi n . 


2nZa 

3 


|T( 0)| 2 



Zam r 

2 


Here, 'F(O) denotes the muon wave function at the origin, 
(r 2 ) is the rms charge radius of the nucleus, and (r 3 )( 2 ) 
is its “Friar moment” El As detailed in Sec. IlII Cl there is 
no contribution from the Friar moment due to a cancel¬ 
lation with part of the inelastic deuteron “polarizability” 
contributions. 

In Eq. Erel and Fnrel contain various moments of 
the nuclear charge distribution, see Ref. [46] , and in par¬ 
ticular the Appendix E therein for analytic expressions 
for some simple model charge distributions. 


The leading order finite size effect, item (rl) in Tab. Ill 


is the first term in Eq. <©• It originates from the one- 
photon exchange with a deuteron form factor insertion 
shown in Fig. [12] and is proportional to the rms charge 
radius of the deuteron, (rjj). 



FIG. 12. Item (rl), the leading nuclear finite size correction 
stems from a one-photon interaction with a deuteron form 
factor insertion, indicated by the thick dot. 

For (r2), the radiative correction a(Za) 5 , we chose 
Martynenko’s value: The equations for the calculation 
of this term are given in m- Borie [221, Tab. 14 uses 
Eq. (10) of Ref. [57] which gives the total radiative cor¬ 
rection of order a(Za) 5 , i.e. the sum of Eqs. (7) and (9) 


4 (E)( 2 ) has been coined “third Zemach moment” by Friar )46] . 
To avoid confusion with the “Zemach radius” rz that appears 
in the finite size effect in the 2S hyperfine splitting (Sec. |Ivc) 
we adopt the term “Friar moment” as recently suggested by 
Karshenboim ESI. 


(f 3 ) (2) + (Za) 2 (F RE L + tn 2 F N REL) 


(9) 


in Ref. [57. Martynenko [23], in contrast, uses Eq. (9) of 
Ref. [52], stating that the additional polarization correc¬ 
tion, Eq. (7) in Ref. 57]. which is included in Eq. (10), 
cancels with a part of the (inelastic) deuteron polariz¬ 
ability contribution. 

The finite size correction to the Lamb shift of order 
(Za) 6 has first been calculated by Friar |46| . see in par¬ 
ticular Appendix E therein. Both Borie [22] (p. 30) and 
Martynenko [24] (Eq. (33)) follow Friar l 46| and evaluate 
this contribution as the sum of two terms which we list 
separately: 

The first one, (r3), has an explicit (r^) dependence, while 
the second one, (r3’), is usually evaluated for an exponen¬ 
tial charge distribution, since a model-independent eval¬ 
uation of this term is prohibitively difficult [22]. Small 
differences between the formulas given by Borie |22j and 
Martynenko [24] result in values for (r3’) of 0.0029 meV 
and 0.0031 meV. For example, the term (7' 2 ) (In(fir )), 
which is part of F RE l in Eq. ([ 9 ]) [^] is attributed to (r3) 
and (r3’) by Martynenko and Borie, respectively. 


TABLE II. The item (r3’), remaining (Za) 6 corrections to 
finite size from Tab. |10[ in meV, evaluated for an exponen¬ 
tial and a uniform charge distribution of the deuteron, using 
the formulas given by Borie [22] p. 30, and Martynenko [21] 
Eq. (33), and the moments from Friar [IB] Tab. V. For de¬ 
tails, see text. Our average, AI?(r3 , ) = 0.0030(6) meV, is 
obtained from the spread of these values. 


Distribution 

after Borie [22 

after Martynenko [21] 

Exponential 

0.00238 

0.00285 

Uniform 

0.00350 

0.00355 


5 See e.g. Ref. 1461 Eq. (43) or Ref. [22] p. 30. 
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TABLE III. Coefficients of the radius-dependent contributions to the Lamb shift. Values are in meV/fm 2 , except for r3’ e . 
KS: Kallen-Sabry, VP: vacuum polarization, SOPT: second-order perturbation theory. 



Contribution 

Martynenko 

Borie 

Karshenboim 


Our choice 




m 

|22| Tab. 14 

130) Tab.Ill 

value 

source 

rl 

Leading fin. size corr., (Za) 4 

-6.07313 

(27) 

-6.0730 

b a 

-6.0732 AE ( ° ] ns 

-6.07310 

±0.00010 

avg. 

r2 

Radiative corr., a(Za) 5 

-0.0009620 

(62), [23] 

—0. 00072[(] 

h 


-0.000962 


M 

r3 

Finite size corr. order (Za) 6 

-0.002128 

(33) 

-0.00212 

b c 


-0.002124 

±0.000004 

avg. 

r4 

Uehling corr. (+KS), a(Za) 4 

-0.01350 

(28) 

-0.0130 

bd 

-0.0132 AEppjg 

-0.01325 

±0.00025 

avg. 

r5 

One-loop VP in SOPT, a(Za) 4 

-0.020487 

(29) 

-0.02062 

be 

-0.0205 AEppjg 

-0.020554 

±0.000067 

avg. 

r6 

Two-loop VP corr., a 2 (Za) 4 

-0.000105 

(30, 31) 




-0.000105 


M 

r7 

Two-loop VP in SOPT, a 2 (Za) A 

-0.000095 

(32) 




-0.000095 


M 

r8 

Corr. to the 2 P \/2 level 



-0.0000606 

fr(2pi/ 2 ) 


+0.00006060 

B 


Sum 

-6.10848 

-6.10952(3 

—6.1069 

-6.11013 

±0.00028 



r3’ 

Remaining order (Za) 6 [meV] 


0.0029 meV (33) 

0.0033 meV 


0.00300 ±0.00060 meV Tab. 

II 



Sum 

-6.10848 r 2 + 0.0029 meV 

-6.10952 r 2 + 0.0033 meV 


-6.11025(28) r d 2 + 0.00300(60) meV 


a This value was published with a wrong sign in m- The term is from Eq. (9) in m 
b This value is obtained in EH, Eq. (10). For further explanations see Sec. 
c The sign is explained in the text. 
d The (r 2 ) coefficient given in Ref. m page 13, neglects the correction to the 2PX/2 level, item (r8). 
e Belongs to r3. Depends on the charge distribution in a non-trivial way, see text. 


IIIB 


We calculate (r3’) from Borie’s and Martynenko’s for¬ 
mula, for both an exponential and a uniform charge dis¬ 
tribution, using the moments given by Friar |46| and ob¬ 
tain the values listed in Tab. [TT] We adopt the average, 
AE(t3') = 0.0030(6) meV. 

The items (r4) and (r5) do not depend on the shape of 
the deuteron charge distribution (221 [35]. The two-loop 
vacuum polarization corrections (r6) and (r7) are only 
given by Martynenko [24]. 

A correction to the 2P ly / 2 level (r8) is given by 
Borie [22] . Item (r8) shifts the 2Pi/2 level “upwards” 


(less bound). This increases the energy difference be¬ 
tween the 2 S and 2 Pi/ 2 levels, which explains the positive 
sign of this contribution in Tab. |III| At the same time, 
this term decreases the fine structure (2P 3 / 2 — 2 P 1 / 2 en¬ 
ergy difference) and is hence listed as item (flO) with a 


negative sign in Tab. VI 


The total radius-dependent contribution to the Lamb 
shift yields 


A^.-dep. = - 6.11025(28) r d meV/fm 2 
+ 0.00300(60) meV. 


(10) 
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C. Nuclear polarizability contributions to the 
Lamb shift 

Historically, the two-photon exchange (TPE) contribu¬ 
tion to the Lamb shift (LS) in muonic atoms has been 
considered the sum of the two parts displayed in Fig. [13] 
(a,b) and (c,d), respectively: 

AE TPE = A-^Priar + AE'inolastic (11) 


Early calculations of the contribution from the 
deuteron polarizability, i.e. the inelastic part AE i I n | lastic 
displayed in Fig. [l3](c,d) include Fukushima etal. 61], 
1.24 meV, Lu and Rosenfelder |S2], 1.45 ± 0.06meV, and 
Leidemann and Rosenfelder [63], 1.500± 0.025 meV. The 
latter value has been used extensively in the literature. 

1. Modern determinations of A1?tpe 


The elastic “Friar moment” contribution, AB^ iar , also 
known as “third Zemach moment” (r 3 )( 2 ) contribution, 
shown in Fig. [T3](a,b) is sensitive to the shape of the nu¬ 
clear charge distribution, beyond the leading (r 2 ) depen¬ 
dence discussed in Sec. |IIIB| This part is traditionally 
parameterized as being proportional to the third power 
of the rms charge radius. The coefficient depends on 
the assumed radial charge distribution. For example, 
Borie gives (r 3 )( 2 ) = 4.0(2) r 3 , where r = i/(r 2 ). For /Ltd, 
the Friar (3rd Zemach) moment contribution amounts to 
~ 0.43 meV [22]. 

The inelastic part, A£’h® lastic , frequently termed “nu¬ 
clear polarizability contribution” is shown in Fig. |13[ c,d). 
It stems from virtual excitations of the nucleus due to the 
exchange of two photons with the muon. The inelastic 
contributions are notoriously the least well-known the¬ 
ory contributions and limit the extraction of the charge 
radius from laser spectroscopy of the Lamb shift. 






FIG. 13. (a)+(b) Elastic AE^ar! and ( c )+(d) inelastic 

AE ia 3 lastic two-photon exchange (TPE) contribution. The 
thick dots in (a) indicate deuteron form factor insertions. The 
blob in (c) and (d) represents all possible excitations of the 
nucleus. The elastic part (a)+(b) is canceled by a part of the 
inelastic (polarizability) contribution of (c)+(d) [55l [58], just 
like for electronic deuterium [SUES]. 


Recently, several works have revisited the TPE contri¬ 
butions to the Lamb shift in /id. Tab. |IV| lists the con¬ 
tributions to A£^p E obtained by Pachucki (2011) 55| . 
Friar (2013) [50] , Carlson etal. (2014) ]M], the 
TRIUMF/Hebrew University group in Hernandez et al. 
(2014) [58], and Pachucki and Wienczek (2015) [65] . 

As it will turn out that the uncertainty in A£pp E is 
by far the largest uncertainty in the determination of 
from the /id data, we next summarize the main features of 
these papers. We identify missing and incorrect terms in 
the original papers. The detailed compilation in Tab. |IV| 
allows us to obtain the reliable average given in Eq. ( [T7[ . 

In his 2011 paper 55], Pachucki calculated the nuclear 
structure corrections to the Lamb shift in muonic deu¬ 
terium using the AV18 potential for the deuteron and ob¬ 
tained AAr^p E = 1.680(16) meV. Moreover, he confirmed 
that for /id, similar to electronic deuterium !sa m, 
the elastic “Friar moment” contribution of order ( Za ) 5 , 
a,b)) is canceled by a part of the inelas- 


AE bLr (Fig- 


13 


tic two-photon (polarizability) contributions, Af?h| 


inelastic 


(Fig.[l3][c,d)). The reason for this cancellation is that the 
deuteron binding energy of 2.2 MeV is small compared to 
the muon mass0 

Pachucki [55] includes both the elastic and inelastic TPE 
contribution of the proton, but not the neutron. For 
the proton, he rescaled the full proton TPE contribu¬ 
tion calculated for muonic hydrogen [66], AE(2S) = 
—0.0369(24) meV, with a reduced mass ratio to correct 
for the larger wave function overlap in /id, 


C = (m^/m^) 3 = 1.1685. (12) 


This gives a value of 0.043(3) meV for our items pl3+pl4. 
Pachucki’s value for the magnetic contribution (p!0) was 


6 For muonic hydrogen, in contrast, the first excited state of the 
nucleus (proton) is the A resonance with an excitation energy 
of 300 MeV. Hence there is no such cancellation between elastic 
and inelastic TPE contributions in /ip. 










12 


found to be wrong by a factor of two in Ref. [58] and was 
corrected in Pachucki’s later work [551 . 

Friar (2013) used the zero-range approximation 
(ZRA) jBOj which allows for a systematic derivation of 
all terms. Friar finds very good agreement with the re¬ 
sults of Pachucki [55] despite the simplicity of the ZRA. 
The cancellation between elastic and inelastic TPE con¬ 
tributions is observed in ZRA, too Q Friar noted that 
a nucleon finite size contribution of 0.029 meV should be 
added that had not been included in Ref. [55] • Friar’s 
value of Ai5^p E = 1.941 ±1% meV seems at first glance 
to be in serious disagreement with Pachucki’s value [55] . 
The difference is however mainly caused by the Coulomb 
distortion (p5+p6) of —0.263 meV which should be in¬ 
cluded in every calculation Including further items 
in Tab. |IV| like the nucleon polarizability contribution 
pl4+pl5, and the nucleon subtraction term pl6 results 
in a “corrected value” of 1.697 meV. Higher order correc¬ 
tions to the dipole contribution (p3 and p4) can account 
for the remaining difference to the other model calcula¬ 
tions and “our avg.”. 

In their 2014 paper [58] the TRIUMF/Hebrew Univer¬ 
sity group performed an independent calculation using 
two parameterizations of the deuteron potential: AV18 
one the one hand, and nucleon-nucleon (NN) forces from 
chiral effective field theory (yEFT) up to order N 3 LO 
and with various cutoffs, on the other. As in their work 
on muonic heliunr-4 [68l [69] they added higher order rel¬ 
ativistic corrections, corrected the magnetic term, and 
added the intrinsic neutron polarizability [70] . They 
also introduced the reduced-mass-dependence in higher- 
order terms, while the earlier Ref. |55j had worked in the 
limit of infinite nuclear mass for all terms. The nuclear 
mass dependence for all terms was then further refined 
in Ref. [65] . 

Ref. [58] observed the cancellation of elastic and inelas¬ 
tic contribution for the muonic deuterium case explic¬ 
itly: The sum of terms <5^/ + = —0.424(3) meV [^] 

in Ref. [58] cancels very nicely with elastic Friar (“3rd 
Zemach”) contribution AU EEar = 0.433(21) meV of the 
deuteron as calculated by Borie, see Ref. [22| p. 7. 
Averaging over their results from AV18 and N 3 LO they 
obtained a value of AA^p E = 1.690(20) meV. The ap¬ 
parently good agreement with Pachucki’s value [55] may 


7 Friar’s paper ESI gives a very good understanding about the 
nature of this cancellation. 

8 The sign convention in Ref. ESI is opposite to the one used here. 


however be accidental as it arises from the cancellation 
of many small differences [551 EDI- Again, adding omitted 
items (pl3 and pl6) results in very good agreement with 
all other sums in Tab. IIVI Note that the 4th column 


in our Tab. IV (“Source 4”) is that value from Tab. 3 
of Ref. [58], columns “N 3 LO-EM” and “N 3 LO-EGM”, 
which deviates most from their AV18 result. This is an 
attempt to be rather conservative when determining “our 
average” following Eq. ([6]). 

tn\ 

According to Bacca [70], their values for S { NS (our item 
pll) should be updated to +0.020meV from the pub¬ 
lished value of +0.015 meV [581. 


The 2015 paper by Pachucki and Wienczek [55] up¬ 
dated Pachucki’s results from 2011 (Ref. [55J), again us¬ 
ing the AV18 potential. Among other things, they in¬ 
cluded the finite size of the nucleons, and the intrinsic 
elastic and inelastic two-photon exchange with individ¬ 
ual nucleons. They corrected their magnetic interaction 
term, and derived the correct mass dependence of the 
TPE correction and its consistent separation with the 
so-called pure recoil correction. Their total TPE con¬ 
tribution of order ( Za) 5 is 1.717 ± 0.020 meV. Item pl6 
must be added to obtain a “corrected value” in very good 
agreement with all other determinations. 


Complementary to the calculations using various 
deuteron potentials [5511581 l60l 165] , Carlson et al. [64] , 
in 2014, determined the TPE contributions with min¬ 
imal model dependence using measured elastic and in¬ 
elastic electron-deuteron scattering data and dispersion 
relations. Their model-independent calculation yields 
2.01(74) meV, confirming the numbers given by [55] [58 , 
I60j albeit with a much larger uncertainty. This uncer¬ 
tainty stems from the uncertainty in the data and can 
be improved significantly when new data from the Mainz 
MAMI and MESA facilities becomes available [B4] , 

The several contributions to their sum can not be 
easily equated with individual items pl...pl6 listed in 
Tab. |IV| so we do not quote their individual contribu¬ 
tions, with one important exception. Carlson et al. [6T 
note that the proton and neutron intrinsic polarizabilities 
of 0.028(2) meV (our items pl4+pl5) should be added to 
the earlier results of Pachucki [55] and Friar [50] , Such a 
correction is already included in the later paper by Her¬ 
nandez etal. [581 . 

On the other hand, the value of Carlson et al. [H] should 
be corrected mm for Coulomb distortion (p5+p6) of 
—0.263 meV. Then the central value becomes 1.748 meV, 
in even better agreement with the (corrected) values from 




13 


nuclear models E3E3ED1E3- Note that this correc¬ 
tion corresponds to 1/3 of their quoted uncertainty and 
may hence look absurd. But the uncertainty quoted in 
Ref. [64] originates almost exclusively from the plain wave 
Born approx. (PWBA) term and may be reduced by at 
least a factor of 4 with new data from a planned exper¬ 
iment in Mainz [64]. The good agreement between the 
corrected central value and all other (corrected) values 
makes one wonder if the uncertainty in Ref. 164] is maybe 
somewhat conservative. 

The “Thomson term” is a recoil correction that has first 
been calculated in Ref. EH. It has received some atten¬ 


tion 03 E3 ED] mi in the discussion of our Tab. pV 


and the conclusion was that this term is indeed cor¬ 
rectly added to the sum of the contributions in the 
dispersion-relation treatment of Ref. [64]. The other 
calculations [55/ ES EDI El] have correctly not included 
such a term, because the cancellation between elastic and 
(part of the) inelastic contributions to the polarizability 
will eliminate this Thomson term (as well as other simi¬ 
lar recoil-like terms) in such a “nuclear Hamiltonian ap¬ 
proach”. All further recoil corrections of order ( Za) 5 to 
the Lamb shift in muonic deuterium are then included in 
the “pure recoil corrections”, item #22 in Tab. p] 


2. Comparison of terms and further corrections 

An earlier version of the present manuscript was sent 
to the authors of Refs. 55, [58l EDI El, 65] and other 
experts in the field. The ensuing insightful discussions 
resolved several discrepancies between the published val¬ 
ues of ALllfp E and revealed that some further corrections 
should be included. 

Table [TV] lists in chronological order the modern deter¬ 
minations of AAlfp E using various nuclear models, and 
scattering data. As usual, we calculate an “average” fol¬ 
lowing our Eq. Q and consider the spread of values in 
the uncertainty. 

Items pi through plO contain the nuclear contribu¬ 
tions, and the various calculations are in good agreement. 

It is satisfying to note that the dominant dipole term, 
item pi, is in very good agreement for the three mod¬ 
els used: AV18, ZRA, and N 3 LO xEFT. We average the 
results using “modern” potentials [13 [65] and take the 
agreement of the ZRA result as an indication that the 
ZRA results for the smaller terms are likely to be accu¬ 
rate on the few /teV level and can hence be used in “our 
average”. 


Items p2..p4 are relativistic corrections to pi. The two 
most recent works [581165] include higher order relativis¬ 
tic corrections so we consider only these works in the 
average. 

There is consensus that the Coulomb distortion con¬ 
tribution (p5+p6) should be included in all calcula¬ 
tions. Adding our average of -0.263 meV to the results 
of Friar [60] and Carlson et al. [64] removes most of the 
discrepancy between all published values. 

Nuclear excitation corrections p7..p9 cancel to some 
degree. Our average includes the results from ZRA [6U 
and the most recent AV18 and N 3 LO models [58] 165] . 

The magnetic contribution plO from Ref. [55] has been 
corrected in the later work [65]. We average over the 
other results. 

Items pll through pl6 are the nucleon contribu¬ 
tions. 

Item pll from the TRIUMF/Hebrew University 
group (6jP s in Ref. [55] ) has been updated to to 
+0.020 meV [ 701] , further improving the agreement with 
Refs. [501 -65]. 

After some discussions, consensus has been reached 
that several nucleon contributions should be included m 
E3E0HZ2]: The elastic Friar (3rd Zemach) term of the 
proton (pl3), the inelastic proton and neutron contribu¬ 
tions (pl4+pl5), and the subtraction terms from both 
the proton and the neutron (pl6) are therefore included 
in our sum. In principle, the elastic Friar term of the 
neutron should be included too, but it is small enough to 
be neglected mmm- 

We follow the suggestion of Birse and McGovern m 
who obtain these values as follows: 

Item pl3, the elastic Friar (3rd Zemach) moment con¬ 
tribution of the proton to the Lamb shift in /id is obtained 
from the values of the elastic and the non-pole Born term 
calculated for muonic hydrogen (/ip) [7T1 . 

Both, the elastic term in /ip, and the non-pole term in /ip, 
have been obtained by Carlson and Vanderhaeghen from 
scattering data using dispersion relations [66!. Their 
value for the elastic term, AE el , amounts to 0.0295 ± 
0.0013meV[/J Their value for the non-pole Born term 
for /ip is —0.0048 meV [/] The sum of these two terms, 
rescaled with ( from Eq. © yields for pl3, the elastic 
Friar moment contribution of the proton to the Lamb 


9 The sign convention in Ref. m is opposite to the one used here. 
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shift in muonic deuterium, the value 

A_£/jjyj ar (p) = 0.0289 ± 0.0015 meV. 


(13) 


The inelastic proton and neutron polarizabilities pl4 
and pl5 have been calculated from deuteron data and 
dispersion relations by Carlson etal. [53] . Their result 
for the sum pl4+pl5 amounts to 


Aplastic(P) + A^i astic (n) = 0.028 ± 0.002 meV (14) 

which is the value we adopt. Hernandez etal. [55, 
used the value 0.027(2) meV from the same Ref. [64] . 
This number is, however, only an estimate using num¬ 
bers rescaled from muonic hydrogen, whereas our choice 


Eq. (14) is calculated from deuteron data, and the value 
in Eq. (14) should be used m- 


Finally, the contribution from the “subtraction term” 
of the nucleon polarizabilities has to be considered, 
too 153 HQ. Birse and McGovern have calculated the 
subtraction term for the inelastic TPE of the proton in 
muonic hydrogen , A E suh = —0.0042 ± 0.0010 meV P*| 
using chiral perturbation theory [73]. This value is in 
good agreement with the value AE subt = —0.0053 ± 
0.0019 meV [^] from Carlson and Vanderhaeghen [66 
which was however obtained from a particular model of 
the proton form factor and an older value of the pro¬ 
ton magnetic polarizability m- For the deuteron, we 
hence adopt the former value, double it assuming that 
the proton and neutron contributions are approximately 


the same El, and rescale with ( from Eq. ( fl2| ) to yield 
pl6 for muonic deuterium 


A£ s L u s b (p) + AE£ b (n) = -0.0098 ± 0.0098meV. (15) 


pLS 


increase conservatively the uncertainty in our average to 
the 0.020 meV obtained by the two most recent model 
calculations [58, 65]. 

The total TPE contribution of order ( Za) 5 to the 
Lamb shift in muonic deuterium □ is hence 

AE^f E (final) = 1.7091 ± 0.0200 meV. (17) 


Further rounding is deferred to Eq. (18). 


The uncertainty of the TPE contribution is by far the 
dominant one, and it limits severely the accuracy of the 
deuteron rms charge radius obtained from laser spec¬ 
troscopy of the Lamb shift in muonic deuterium. 


Here we have assigned a 100% uncertainty. 


3. Our choice 

Summing all values in Tab. m and adding the uncer¬ 
tainties from (the spreads of) our averaging in quadra¬ 
ture, gives 

AE^| e (simple) = 1.7091 ± 0.0146 meV. (16) 

This uncertainty is smaller than the published uncertain¬ 
ties in all original papers [551 [55] [501[65]. Hence we 


11 Note that non-perturbative Coulomb corrections of higher order 
in ( Za ) have been accounted here, whereas the pure recoil part 
of the TPE has been separately given in Tab. IH #22. 


10 The sign convention in Ref. [731 is opposite to the one used here. 














TABLE IV. Deuteron structure contributionsto the Lamb shift in muonic deuterium. Values are in meV. For source 4, the N 3 LO ^ calculation by Hernandez etal. |58| 
we use their value from the rightmost two columns of their Tab. 3 that differs most from their “AV18” value. Items with a diamond ^ are corrected from the published 
values, see footnotes. 
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D. Total Lamb shift in muonic deuterium 


Collecting the radius-independent (mostly) QED con¬ 
tributions listed in Tab. [I] and summarized in Eq. ([8]), 


the radius-dependent contributions listed in Tab. Ill and 
summarized in Eq. (101, and the complete two-photon 
(polarizability) contribution from Eq. (17), we 

obtain for the 2S — 2P energy difference in muonic deu¬ 
terium 


AE(2S-2P 1/2 ) = 228.77356(75) meV 

+ 0.00300(60) meV - 6.11025(28) r 2 A meV/fm 2 
+ 1.70910(2000) meV 

= 230.486(20) meV - 6.1103(3) r 2 A meV/fm 2 


(18) 


where in the last step we have rounded the values to 
reasonable accuracies. 

One should note that the uncertainty of 0.020 meV 
from the nuclear structure corrections AE^| e , Eq. (l7| , 
is about 30 times larger than the combined uncertainty of 
all radius-independent terms summarized in Tab. [TJ and 
15 times larger than the uncertainty in the (r 2 ) coefficient 
(which amounts to 0.0013meV). A further improvement 
of the nuclear structure corrections in light muonic atoms 
is therefore desirable. 


IV. 2S HYPERFINE SPLITTING 

A. Fermi and Breit contributions 


The interaction between the magnetic moment of the 
nucleus with the magnetic field induced by the lepton 
gives rise to shifts and splittings of the energy levels 
termed hyperfine effects. In classical electrodynamics, 
the interaction between the magnetic moments fi A and 
of deuteron and muon, respectively, is described by 

m 


rrclassical 

^HFS 


= -gMd • P^(r) 


(19) 


where S(r) is the delta-function in coordinate space. A 


similar Hamiltonian to the one in Eq. (19) can be derived 


in quantum field theory from the one-photon exchange 
diagram. Using the Coulomb wave function, this gives 
rise in first-order perturbation theory to an energy shift 


for muonic deuterium nS-states of 


r-i 4 (Za) 4 ml 1 

E rfs (F) = x-g-(1 + k)( 1 +a M )- 


— g AUpermi 


F(F + 1) — 


11 


F(F + 1)-^ 


( 20 ) 


where AEy errn \ is the Fermi splitting, m A is the deuteron 
mass, F is the total angular momentum, k and a^ are 
the deuteron and muon anomalous magnetic moments, 
respectively. 

The Fermi splitting 


with 


2 {Za) A m^. 

AApermi = -;-(1 + At)(l + Op) 

rv i m^m A 

= 2^d( 1 + °/i) 


h b = 1 ^M. a + k) 


( 21 ) 


?>n?m ll m A 


( 22 ) 


is the main contribution to the HFS, (hi) in Tab.[V| The 
value Borie gives on p. 19 of Ref. 


is 


AU Eermi = 6.14298 meV. 


(23) 


It already includes the correction AE^amm (h4) due to 
the muon anomalous magnetic moment (/iAMM). How¬ 
ever, Borie’s value AU Eermi = 6.14298 meV is not cor¬ 
rect. Since AEy ernA depends only on fundamental con¬ 


stants @], we have recalculated it following Eq. (21), and 
obtain a value of 


AA Fermi = 6.14308 meV (24) 

which differs from Borie’s, but coincides with Marty¬ 
nenko’s value 
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A E¥ e - = 6.1431 meV 

( 25 ) 

= 6.1359 meVhi+0.0072 meV h4 . 

Here, the two terms in the second line are, respectively, 
the Fermi splitting excluding the contribution of the 
anomalous magnetic moment of the muon, (hi), and the 
/iAMM correction, (h4), which Martynenko calculates 
separately. 

The Breit term AEpreit (h2) corrects for relativistic 
and binding effects accounted for in the Dirac-Coulomb 
wave function but excluded in the Schrodinger wave func¬ 
tion. Both, Martynenko and Borie, calculated the Breit 
correction term to be 


There is no considerable contribution from weak interac¬ 
tion [73] . 

C. Zemach radius 

The Bohr-Weisskopf effect m is the main finite size 
correction to the 2S hyperfine splitting. It is also called 
the Zemach term US], A ^Zemach> and is listed as item 
(h20) in our summary. The Zemach term is usually pa¬ 
rameterized as m 

A^Zemach = -AE Fermi 2 (Zot)m r T Z (29) 
using the so-called Zemach radius of the nucleus 


A^sreit = 0.0007 meV. (26) 


B. Vacuum polarization (VP) and self-energy (SE) 
contributions 

On p. 21 of Ref. [22] , Borie provides the values for eypi 
(h8) and eyp 2 (h5). The values we used are the ones for 
a point-like nucleus. In Sec. |IV C'| we give a correction 
due to the finite size, which is given by (h25, h26). The 
corresponding terms from Martynenko and Borie agree. 

(h7) is neglected by Martynenko as pointed out on p. 21 
in [22]. The origin of the difference between Martynenko 
and Borie in (h9) is not clear. In this case we take the 
average. (h9b) is a correction in third order perturbation 
theory which is only given by Martynenko. 

The pVP contribution AE^yp (hl2) is given by Mar¬ 
tynenko as 

AE^yp = 0.0002 meV. (27) 

Borie included this contribution in the vertex term 
AE vert ex (hl3) as pointed out on p. 21 in [22]. As we 
do not use Borie’s vertex term, we have to consider the 
pVP term from Martynenko. 

Martynenko gives a value for a term called radiative 
nuclear finite size correction (hl7b). It is composed of 
four terms, pSE with nuclear structure , jellyfish correc¬ 
tion and two vertex correction terms, so it should also 
include (hl3). Their sum yields —0.0005 meV. We think 
that (hl4) is an additional term, only calculated by Borie. 
So we add this to ’our choice’. 

The hadron VP AEhyp (hl8) results equal for both, 
Martynenko and Borie: 

(28) 



dV r'p E (r)p M (r - r') 


=- [ (G E (q 2 )G M (q 2 ) - l) ■ 

7T Jo q 2 


(30) 


This convolution of charge p E (r) and magnetization 
Pm(t) distribution comes from the fact, that the finite 
charge distribution alters the muon’s wave function at 
the origin. 

Diagrammatically, the Zemach correction (h20) to the 
HFS, AE£f m s ach , is the elastic part of the two-photon ex¬ 
change contribution to the 2S HFS, just like the Friar 
correction to the Lamb shift Fig. [l3](a,b) [IS] . 

A recoil correction to the elastic TPE, (h23), is consid¬ 
ered here, too. It is somewhat parallel to the item ff22 
of the Lamb shift, termed rel. RC ( Za) 5 listed in Tab.JI] 
The inelastic part of the TPE correction to the 2S HFS, 

A^inellstiO is t0 P ic of the next Sec - 

Borie gives a value of ez emac h = —0.007398 fm -1 r z 
(Ref. []22] , P-22). Her Zemach contribution is hence 
(Eq. @, Ref. [22], P- 23 top) 

AA^nach (Borie) = -0.04545 meV/fm. (31) 


IV D 


Note that the coefficient 

- 0.04545 = | p D (1 + a M ) e Zema ch (32) 

does explicitly include the factor (1 + a^). 

Using r z = 2.593 ± 0.016 fm from Ref. [77] Borie’s value 
for the Zemach contribution to the 2S HFS amounts to 


A£fem ach (Borie) = -0.11782 ± 0.00074 meV. (33) 

Borie mentions that nuclear recoil corrections (our item 
(h23)) are important, but have not been included 
(Ref [22] , P- 22, bottom). 


AE^to' = 0.0002 meV. 
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Martynenko calculates the nuclear structure correction 
a 5 from the deuteron electromagnetic current that in¬ 
volves the form factors F\, F 2 , and F 3 which can be re¬ 
lated to the measured charge, magnetic and quadrupole 
form factors of the deuteron, Ge(q 2 ), Gm(q 2 ) and 
Gq(< 7 2 ), see Eq. (37) in Ref. [15]. Using the parameter¬ 
ization of Ge, Gm and Gq from Ref. [78], Martynenko 
obtains a value of — 0.1163±0.0010 meV (Ref. [15], Tab. 1 


item #7, Eq. (46)). 

This value is the sum of the Zemach term, item (h20), 
as calculated by Borie, but includes recoil corrections to 
the finite size effect, (h23) (331 ■ The separation of these 
two contributions is not unique [33], but if one adopts 
the canonical definition of the Zemach radius in terms of 


the form factors Ge and Gm as given in Eq. (30), one 
can separate Martynenko’s sum into 


A ^emach+RC (Martynenko) 


- 0.1163 meV ± 0.0010 meV 

(—0.1178 meVh 20 + 0.00 1 5 meVh23) ± 0.0010 meV. 


(34) 


r 


Here the item (h20) was calculated using the Zemach ra¬ 
dius rz = 2.5959 fm [33], obtained by numerical integra¬ 
tion of the parameterization of the deuteron form factors 
from Ref. [73]. This Zemach radius is in excellent agree¬ 
ment with the value rz = 2.593 ± 0.016 fm from Ref. 
used by Borie [22] . A small difference arises from Marty¬ 
nenko’s observation that the factor (1 + a M ) in Eq. (21) 


should not be included for the 2 y amplitudes with point 


vertices 155] . 

Rewriting the — 0.1178 meV of (h20) as 

AZemach (Martynenko) = —0.0453934 r z meV/fm 

(35) 

makes the dependence on the Zemach radius explicit. 
Putting the nuclear recoil corrections back in, the com¬ 
bined Zemach (h20) and recoil (h23) corrections evalu¬ 


ated by Martynenko (Eq. (34)) become 


A£ fe^ach+Rc = -0.0453934 rz meV/fm h20 + (0.0015 ± 0.0007) meV h23 (36) 


which we adopt. The total uncertainty of 0.0010 meV 
given by Martynenko is then the sum of the uncertainty 
in the Zemach radius Srz = 0.016 fm, corresponding to 
0.0007 meV, and the uncertainty of (h23) given above. 

D. Nuclear polarizability contributions to the 2S 
HFS 

The polarizability contribution to the 2S hyperfine 
splitting in muonic deuterium, AUfjfpf, has only recently 
been calculated for the first time by the group of 
Martynenko M- They obtain the polarizability term in 


two parts: 

First, the deuteron polarizability contribution 
AE™! (deuteron) (h 22 a) is obtained from the ana¬ 
lytic expressions derived in zero range approximation for 
electronic deuterium by Khriplovich and Milstein m- 
This part takes into account the virtual excitation of a 
deuteron made from point nucleons. 

Second, the much smaller internal deuteron polarizabil¬ 
ity contribution Ai?i nt . d-poi. (h22b), which accounts for 
the excitation of the individual nucleons (proton and neu¬ 
tron) inside the deuteron. This part is estimated based 
on the results for muonic hydrogen [80] . 

Summing these two up yields 


A ^tot. d-poi. = 0.2121(42) meV + 0.0105(25) meV = 0.2226(49) meV (37) 


with a generous uncertainty that accounts also for the fact that the original derivation ;79| was for electronic, 
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and not for muonic deuterium. 


Eq. (37) is also the value quoted by Borie. As pointed 
out by Borie, Ref. [22] p. 22, it is not clear whether the 
’elastic’ contribution of the two-photon exchange dia¬ 
grams is taken into account. 

As for the Lamb shift, the polarizability term for the 
2S HFS is the one with by far the largest uncertainty. 


E. Further corrections to the 2S HFS 

Several further corrections are considered by either 
Martynenko or Borie (h24, h25, h26, h27, and h27b): 
Martynenko calculates a mixed term which includes eVP 
and a nuclear structure correction (h24) 

A£ eVP+ nucl . str u Ct . = 0.0019 meV (38) 

as well as two further nuclear structure contributions 
(h27, h27b): 

A< cl , tr , on , = 0.0008 meV (39) 

A< ucl . str . SO PT = -0.0069 meV (40) 


Borie gives finite size corrections (h25, h26) to the eVP 
terms eypi and eyp 2 ) both contributing -0.00068meV. 
These are obtained by calculating the difference of eypi 
and eyp 2 with a point size nucleus compared to the ones 
when considering a finite size ([25], p. 21). It is not yet 
clear whether the nuclear structure contributions from 
Martynenko are complementary to the ones of Borie. We 
should remark that item (h27b) is quite big. It doesn’t 
seem to be included in Borie’s calculations and is the 
main reason for the difference between Borie 1221 and 
Martynenko 1151 . 

For now we refrain from assigning a large uncertainty 
to this item h27b, but an independent calculation, or 
at least an estimate of it’s accuracy, would certainly be 
helpful. 


F. Total 2S hyperflne splitting 


Hence, collecting all terms, but separating out the 


deuteron polarizability correction Eq. (37) as it is the 


dominant source of uncertainty, we can write the total 
2S HFS in muonic deuterium as 


A£$fs = 6.17415(73) meV + 0.22260(490) meV - 0.04539 meV 
= 6.39675(494) meV - 0.04539 meV. 


(41) 


The large uncertainty in the polarizability corrections to 
the 2S HFS will prevent a determination of the deuteron 
Zemach radius from the measured transitions in muonic 
deuterium m- An improved calculation of the polariz¬ 
ability terms is therefore highly desirable. 

Using the Zemach radius rz = (2.593 ± 0.016) fm [77] 
we get: 


to be compared to the muonic deuterium measure¬ 
ment 1101 . Alternatively, one can use the measurement 
and the Zemach radius to accurately determine the po¬ 
larizability contributions. Such a number may serve as a 
benchmark for accurate lattice calculations. 


AF'hfs = 6.27905(495) meV 


(42) 
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TABLE V. All contributions to the 2S hyperfine splitting (HFS) in muonic deuterium. The item numbers hi in the first 
column follow the entries in Tab. 3 of Ref. [3]- For Martynenko, numbers #1 to # 15 refer to rows in Tab. I of Ref. m, whereas 
numbers in parentheses refer to equations therein. Borie |22l gives the values as coefficients to be multiplied with the sum 
of (hl+h4). We list the resulting values in meV. AMM: anomalous magnetic moment, PT: perturbation theory, VP: vacuum 
polarization, SOPT: second order perturbation theory, TOPT: third order perturbation theory. 


All values are in meV (meV/fm for h20). 



Contribution 

Martynenko [15] 

Borie [22] 


Our choice 


hi 

Fermi splitting, ( Za) 4 

6.1359 

#L ( 6 ) 





h4 

pAMM corr., a(Za ) 4 

0.0072 

#2, (7) 




Eq. J24| 

sum 

(hl±h4) 

6.1431 



6.14298 

p. 19 

6.14308 

h 2 

Breit corr., (Za) 6 

0.0007 

#3, ( 8 ) 

0.00069 

p. 19 

0.00069 

B 

h5 

eVP in 2nd-order PT, a(Za) 5 (£yP 2 ) 

0.0207 

#4, (23) 

0.02070 

p. 21 

0.0207 

M 

h7 

Two-loop corr. to Fermi-energy (eyp 2 ) 

neglected 



0.00016 

p. 21 

0.00016 

B 

h 8 

One-loop eVP in I 7 int., a(Za) 4 (eypi) 

0.0134 

#4, ( 12 ) 

0.01339 

p. 21 

0.0134 

M 

h9 

Two-loop eVP in I 7 int., a 2 (Za) 4 (eypi) 

0.0005 

#5, (16), (29-32) 

0.00010 

p. 21 

0.0003 ±0.0002 

avg. 

h9b 

VP contr. in TOPT 

0.00004 

# 6 , (33) 



0.00004 

M 

hl 2 

fjiVP (sim. to eyp) 

0.0002 

#9, (48) 

incl. in hl3 

p. 21 

0.0002 

M 

hl3 

Vertex, a(Za) 5 

inch in hl7b 



-0.00059 

p. 19 

incl. in hl7b 


hl4 

Higher order corr. of (hl3), part with ln(a) 




-0.00004 

p. 19 

-0.00004 

B 

111 7b 

Radiative nucl. fin. size corr., a(Za) 5 

-0.0005 

#13, (71-74) 



-0.0005 

M 

111 8 

Hadron VP, a 6 

0.0002 

#10, (50) 

0.00016 

p. 19 

0.00016 

B 

hl9 

Weak interact, contr. 

0 

p. 10 


0 

p. 21 

0 

mi 

h 20 

Fin. size (Zemach) corr. to A-Epermij (Za) 5 

-0.04539 r z 

ES 


-0.04545 r-43 

p. 22 

0.04539 r z 

M 

h23 

Recoil corr. to fin. size 

0.0015 ±0.0007 

m 




0.0015 ±0.0007 

M 

sum 

(h20+h23) 

-0.1163 ±0.0010 

#7, (46) 




M 

h 22 a 

Deuteron polarizability, (Za) 5 

0.2121 ±0.0042 

#14 using [79] 



0.2121 ±0.0042 

M 

h 22 b 

Deuteron internal polarizability, (Za) 5 

0.0105 ±0.0025 

#15 using (80 



0.0105 ±0.0025 

M 

sum 

(h 22 a+h 22 b) 

0.2226 ±0.0049 



0.2226 ±0.0049 

p. 22 



h24 

eVP + nucl. struct, corr., a 6 

0.0019 ±0.00001 # 8 , (47) 



0.0019 

M 

h25 

eVP corr. to fin. size (sim. to eyp 2 ) 




-0.00068 

p. 21 

- 0.00060 

B 

h26 

eVP corr. to fin. size (sim. to eypi) 




-0.00068 

p. 21 

-0.00068 

B 

h27 

Nucl. struct, corr., a(Za) 5 

0.0008 

#11, (55) 



0.0008 

M 

li27b 

Nucl. struct, in SOPT 

-0.0069 

# 12 , 

(59) 



-0.0069 

M 


Sum 

6.39824 ±0.00494 



6.39880 ±0.00490 


6.39675± 0.00494 




-0.04539 r z 



-0.04545 r z 


-0.04539 r z 



a The published value for the sum of items h20+h23 is —0.1163 ± 0.0010 meV 1151 . For the separation into items h20 and h23 see text. 
b Calculated from Eq. ( |2l| , including the factor (1 + a M ). According to Martynenko, this factor should be omitted in the 2c amplitudes 
with point vertices. 

c Difference of two terms in Borie m- See text. 
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TABLE VI. Contributions to the 2P fine structure. The items (f7a), (f7d), and (f7e) originate from the same graphs as 
the Lamb shift items #11, #12, and #30*, respectively. VP: vacuum polarization, AMM: anomalous magnetic moment, KS: 
Kallen-Sabry. 

All values are in meV. 



Contribution 

Martynenko 

[241 Tab. 2 

Borie 

M 

Karshenboim 

Our choice 

ft 

Dirac 


8.86430 




f2 

Recoil 


-0.02521 




f3 

Contrib. of order ( Zol) 4 

8.83848 





f4 

Contrib. of order ( Za) 6 and ( Za) Q m\/m 2 

0.00030 





sum 

(fl+f2) or (f3+f4) 

8.83878 

8.83909 


8.83894 ± 0.00016 avg. 

f5 

eVP correction (Uehling), a(Za) 4 

0.00575 

0.00575 

0.0057361 [30] Tab.IV 

0.0057361 K 

f6 

2nd order eVP corr. (KS), a 2 (Za) 4 

0.00005 

0.00005 

0.0000501 [37] Tab.IX “eVP2” 

0.0000501 K 

f7a 

a 2 (Za) 4 m, like #11 



0.0000127 [37] Tab.IX (a) 

0.0000127 K 

f7d 

a 2 (Za) 4 m, like #12 



0.0000991 [33 Tab.IX (d) 

0.0000991 K 

f7e 

a 2 (Za) 4 m, like #30* 



0.0000012 [37] Tab.IX (e) 

0.0000012 K 

f8 

AMM (second order) 


0.01949 




f9 

AMM (higher orders) 


0.00007 




sum 

Total AMM (f8+f9) 

0.01957 

0.01956 


0.019565 ±0.000005 avg. 

flO 

Finite size 

-0.00028 

-0.00027 


-0.000274 



Sum 

8.86387 

8.86419 


8.86412 ±0.00016 


a This is item (i'8), evaluated for the deuteron radius from Eq. J51, see text. 


V. 2P FINE STRUCTURE 

The contributions to the 2P 3 / 2 —2Pi/ 2 fine structure 
splitting in muonic deuterium are displayed in Tab. m 

The main contributions to the fine structure have 
only been calculated by Borie [22] and Martynenko’s 
group (23j [24] • For the latter, we refer to the more recent 
paper Ref. [24] . The values agree in both papers. 

As always, Borie starts from the Dirac equation, which 
has to be corrected for recoil effects. This sum of entries 
(fl)+(f2) has to be compared to Martynenko’s leading 
term (f3), corrected for relativistic effects (f4) which are 
automatically included in the Dirac equation. The result 
of both approaches agree reasonably well and we adopt 
the average 8.83894 ± 0.00016 meV. 

The relativistic recoil correction of order a(Za) 4 , (f5), 
has been calculated including all recoil corrections of or¬ 
der to/M by Karshenboim [30] . This value thus super¬ 
sedes [35] the values obtained by Borie [22] and Marty¬ 
nenko [23l [24l . 

The Kallen-Sabry term, item (f6), agrees nicely among 
all authors. Karshenboim et al. have evaluated some 
higher order a 2 (Za) 4 m contributions with great accu¬ 
racy m which we list as (f7a), (f7d), and (f7e). These 
terms originate from the same graphs as the Lamb shift 


items #11, #12, and #30*, shown in Figs. [6] [7] and [9] 
respectively. 

Contributions from the anomalous magnetic moment 
of the muon in second (f8) and higher (f9) orders have 
been calculated by Borie and Martynenko et al., and their 
sums agree. 

A finite size correction to the 2 P 1 / 2 state (flO) is in¬ 
cluded, too. This term is the same term as item (r8) of 
the (r 2 )-dependent contributions to the Lamb shift, but 
with the opposite sign (see discussion in Sec. |IHB] ). We 
evaluate (r8) for the deuteron radius from Eq. © and 
obtain for our choice of the item (flO) —0.000606 (r 2 ) = 
—0.000274 ± 6 • 10 -8 meV. 

Summing up we obtain as our choice for the 
2 P 3 / 2 - 2 P 1 /2 fine structure splitting in muonic deuterium 

A£ fs (2P 3/2 - 2P 1/2 ) = 8.86412(16) meV. (43) 

VI. 2P LEVELS 

The various 2P levels displayed in Fig. |T] are separated 
by the 2P fine structure treated in Sec. [Vj and further 
split by the 2P hyperfine splitting caused by the magnetic 
hyperfine interaction and the electric quadrupole inter¬ 
action. The Breit-Pauli Hamiltonian can be displayed in 
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matrix form as a sum of the magnetic HFS matrix and the quadrupole interaction matrix: 



/-l.381777 

0 

-0.126405 

0 

0 > 



0 


0.690889 

0 

-0.199864 

0 


MBreit—Pauli — 

-0.126405 

0 

8.161148 

0 

0 



0 


-0.199864 

0 

8.582931 

0 



V 0 


0 

0 

0 

9.285903/ 

magnetic HFS 


/ 0 


0 

0.613872 

0 

° \ 



0 


0 

0 

-0.194123 

0 


+ 

0.613872 

0 

0.434073 

0 

0 



0 


-0.194123 

0 

-0.347258 

0 



V 0 


0 

0 

0 

0.086815/ 

quadrupole int 



/ 

9 p^ = l / 2 
Z - r l /2 

9 p^ 1 — 3/2 
Zr i /2 

9 p^ = l / 2 
^3/2 

9 p-^ 1 —4/2 
Zr 3/2 

9 pF=5 /2 
Zr 3/2 

0 

op^=l /2 

jLr \/1 

-1.381777 

0 

0.487467 

0 

9 pF—3 /2 
Zi l /2 


0 

0.690889 

0 

-0.393988 

0 

= npF=l/2 

Zr 3/2 


0.487467 

0 

8.595220 

0 

0 

9 pF = 3 /2 
Zr 3/2 


0 

-0.393988 

0 

8.235672 

0 

9 pF=5 /2 
Z - r 3/2 

V 

0 

0 

0 

0 

9.372718^ 


It attains off-diagonal elements from mixing of levels 
with the same total angular momentum F, but differ- 


TABLE VII. Elements of the transition matrix (Eq. (441) from 
Ref- [22] . See also Ref. [81] , The various quantities are ex¬ 
plained in the text. 


j j' Energy 
magnetic HFS 

1/2 1/2 (/3/j/6)(2 + Xd + a M )[—<5 f,i/2 + 1/2 <5f,3/2] 
3/2 3/2 AEf a + (^/4)(4 + 5 x d - a^x 


[—l/6<5i?i/2 — 1/15(5^3/2 + 1/10 5 Fi 5/2] 
3/2 1/2 (j3' n /48)(l + 2x d -a fl )[-V2S Fil/2 - V55 Ft3/2 ] 


quadrupole interaction 

1/2 1/2 0 

3/2 3/2 tQ [5 f, 1/2 — 4/5 Jp,3/2 + 1/5 (Si? 5/2] 
3/2 1/2 eq [V25 Fi 1/2 -1/V55 Fi3/2 \ 


ent total muon angular momentum j di m m mi , as 
shown in Tab. |VII| Note that the diagonal terms of the 
quadrupole interaction lead to a change in the order of 
the 2P 3 /2 levels (see also Fig. [l]). 

We follow Borie’s treatment [22], see also Pachucki [15] 
and Jentschura [3Bj, but use our value for the 2P fine 
structure AEf s (2P 3 / 2 — 2 Pi/ 2 ) = 8.86412(16) meV from 
Sec. [VJ Eq. ( |43| , as well as a more recent value of the 
deuteron quadrupole moment 


Q = 0.285783 (30) fm 2 


(45) 


from Ref. [52], or, equivalently, Q = 7.33945(77) ■ 
10” 24 /meV 2 , using he = 197.3269718(44) MeV fm g]. 
The numerical values of the quantities used in Tab. |VIl1 
are given in Tab. VIII In brief, /3p is defined in Eq. (221. 
For the 2P levels, 


P'd = /3d(1 + £2p) (46) 

has to be used which contains the Uehling correction re¬ 
quired for levels with l > 0 (see p. 25 and Eq. (12) in 
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TABLE VIII. Input parameters for the transition matrix. We 
recalculated Borie’s values, but use our fine structure (see 
Sec •0 and an updated value of the quadrupole moment Q 
[82l for tQ , see Eq. (481. 



Borie [22] 

our value 

Pd 

4.0906 meV 

4.0906259 meV 

Pd 

4.0922 meV 

4.0922253 meV 


0.43439 meV 

0.434073(46) meV 

x d 

0.0248 

0.0247889 


0.00116592 

AE{ S 

8.86419 meV 

8.86412(16) meV 


Ref. [22])- For muonic deuterium, 


e 2P = 0.000391. (47) 


The quadrupole moment of the deuteron enters in the 
hyperfine splitting via the quadrupole interaction, see 
Borie [22], pp. 24 and 25. 


e Q 


aQ ( aZm r ) 3 
~2 24 


(1 + ^2 p), 


(48) 


where Q the quadrupole moment of the deuteron and e 2 p 
is given in Eq. (471. 


Diagonalizing the matrix Eq. (44) results in shifts of 
the 2P(F = 1/2) and 2P(F = 3/2) levels by 


Ai /2 = 0.02376 meV and 
A 3 / 2 = 0.02052 meV, 


(49) 


respectively, as displayed in Fig. [I] The resulting energies 
of the various 2P sublevels are summarized in Tab. IIXI 


TABLE IX. 2P levels from fine- and hyperfine splitting. All 
values are in meV relative to the 2P!/ 2 level. The fine struc¬ 
ture (2P 3 /2—2 Pj/2 energy splitting) is our value AE{ S = 
8.86412(16) meV from Eq. (431. Uncertainties arise from the 
quadrupole moment Q in Eq. (451 and AEc s . 


Borie [2!| Our value 


9 p^ 1 —1/2 

Zr l/2 

-1.4056 

-1.40554(1) 

9 p F= 3/2 
Z ' r i/2 

0.6703 

0.67037(1) 

9 TD ^ = 1/2 

Zr 3/2 

8.6194 

8.61898(17) 

9 p F= 3/2 
Zr 3/2 

8.2560 

8.25619(16) 

op F= 5/2 
Zr 3/2 

9.3728 

9.37272(16) 


VII. SUMMARY 


In summary, we have compiled all known contributions 
to the Lamb shift, the 2P fine structure, and the 2S and 
2P hyperfine splittings, from QED and nuclear structure 
contributions. 

For the Lamb shift, the QED contributions in Tab. |T] 
show good agreement between the four (groups of) au¬ 
thors. A problem with our item #2 from Ref. ;25] was 
identified and resolved by the auhors. Ultimately, the 
uncertainty of these “pure QED” terms in Tab. [j] is suf¬ 
ficiently good. 

For the radius-dependent terms in Tab. |III| we find good 
agreement between the authors, too. Some terms have 
however been calculated by only one group. We re¬ 
calculate a small term (r3’) to verify that the model- 
dependence imposed by this contribution is sufficiently 
small. 

The main limitation for the Lamb shift, and hence the 
deuteron charge radius to be extracted from the upcom¬ 
ing data, originates from the two-photon exchange con¬ 
tribution to the Lamb shift in /id. Here, a superficial in¬ 
spection of the six modern values published in five papers 
since 2011 ESI ESIMI EH vary between 1.68 meV and 

2.01 meV, suggesting an uncertainty as large as 0.3 meV. 
The term-by-term comparison of the individual contri¬ 
butions in Tab. m revealed that the agreement is in 
fact much better. Fruitful discussions with the authors 
of these papers and other experts in the held revealed 
missing terms and resulted in updated values of some 
individual terms. It is very reassuring that vastly differ¬ 
ent approaches give results in excellent agreement, when 
corrected for missing terms: zero-range approximation, 
modern nuclear models like AV18 (from two groups of au¬ 
thors) and \EFT-inspired NN-forces up to N 3 LO order, 
and dispersion relations using electron-deuteron scatter¬ 
ing data. Our average, 1.709 ± 0.020 meV is a reliable 
prediction for the deuteron polarizability contribution to 
the Lamb shift in /id. 

For the 2S-HFS, several nuclear structure contribu¬ 
tions have so far only been calculated by one group m-- 
These are items (h22a), (h22b), and (li27b) in Tab. 0 
which are rather large, and their uncertainties dominate 
the theoretical uncertainty for the 2S-HFS. This uncer¬ 
tainty will prevent us from obtaining a meaningful value 
of the Zernach radius of the deuteron from the measure¬ 
ment of the 2S-HFS in /id. An improved calculation of 
these items is therefore desirable. 
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For the 2P fine- and hyperfine splittings we collect all 
terms from the various authors, recalculate the matrix el¬ 
ements of the Breit-Pauli Hamiltonian with updated val¬ 
ues of the 2P fine structure and the deuteron quadrupole 
moment. Diagonalizing this matrix Eq. (441 we obtain 
the 2P level energies, and their uncertainties. 
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Appendix A: The Darwin-Foldy term 

The Darwin-Foldy term, which is part of the Barker- 
Glover corrections (our item #17 in Tab. |T| has histori¬ 
cally been subject of different definitions. 

Pachucki and Karshenboim [49] argue that the DF 
term originates from the Zitterbewegung of the nucleus 
and is hence absent for a spin-1 nucleus such as the 
deuteron (as well as for the spin-0 4 He nucleus). 

Khriplovich, Milstein and Sen’kov m argue that the 
DF term must be made a part of the rms charge radius 
to be consistent with electron scattering. In this case, 
the DF term is not absent for spin-1 nuclei such as the 
deuteron. 


Friar, Martorell and Sprung m have emphasized that 
the DF-term can be alternatively considered as part of 
a recoil correction of order 1/M 1 2 , or as the energy shift 
due to a part of the mean-square radius of the nuclear 
charge distribution. They advocate the second choice but 
admit that the first choice has to be used for the proton 
because ”it is unfortunately far too late to change these 
conventions for the hydrogen atom”. They recommend, 
however, to not extend the hydrogen atom conventions 
to other nuclei. 

Jentschura has discussed the situation in some 
breadth [52] and concluded that the DF term should in¬ 
deed be considered a contribution to the atomic energy 
levels due to the nuclear Zitterbewegung, supporting the 
“atomic physics” convention of Ref. 4S! . The DF term 
is hence absent for the deuteron. 

This “atomic physics” convention, in which the DF 
term is not a part of the rms charge radius, but rather 
a recoil correction of order (Za) A m 3 /M 2 to the energy 
levels, is the convention used in CODATA-2010 [3], see 
Eq. (26) and (27) therein. It is also the convention used in 
the most recent measurement of the H-D isotope shift m 
EH which is the origin of the difference of the squared rms 
radii of the deuteron and the proton given in Eq. 
Moreover, it is the convention used for the proton radius 
in muonic hydrogen [fl . 

Therefore, to be able to directly compare the numer¬ 
ical values of the proton and deuteron rms charge radii 
obtained in electronic and muonic atoms, one must follow 
the “atomic physics” convention gHUilllGS], which is 
what we do. 
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